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Abstract 

(N 

A study on filtering aspects of airborne wind energy generators is presented. This class of renewable energy systems 
aims to convert the aerodynamic forces generated by tethered wings, flying in closed paths transverse to the wind flow, 
£\J ■ into electricity. The accurate reconstruction of the wing's position, speed and heading is of fundamental importance 

for the automatic control of this kind of systems. The difficulty of the estimation problem comes from the nonlinear 
dynamics, wide speed range, large accelerations and fast changes of direction that the wing experiences during operation. 
It is shown that the overall nonlinear system has a specific structure allowing its partitioning into sub-systems, hence 
■ leading to a series of simpler filtering problems. Different sensor setups are then considered, and the related sensor fusion 

algorithms are presented. The results of experimental tests carried out with a small-scale prototype are discussed. The 
CNj ■ designed filtering algorithms rely purely on kinematic laws, hence they are independent from features like wing area, 

aerodynamic efficiency, mass, etc. Therefore, the presented results are representative also of systems with larger size and 
, different wing design, including systems with one and two tethers and/or rigid wings. 

r/5 '. 1 Introduction 

High-altitude wind is a vast, still untapped renewable source of energy that has received an increasing attention in the last 
decade, from both industry and academia, with the ongoing development of a series of technologies that fall under the 
umbrella name of airborne wind energy. While many concepts of airborne wind energy generators were present already 
in patents and publications in the late 1970s ll2Tl[T6l[T9l , it is only in recent years that an increasing number of research 
groups and companies started to develop operating prototypes to convert the energy of high-altitude wind, blowing up to 
1000 m from the ground, into electricity (see e.g. [23] [T] [2l [3] [71 [121 El E4] [3 [25]| ) . The activities carried out worldwide 
in the last six years allowed to define and assess, through theoretical, numerical and experimental research, some common 
grounds of airborne wind energy (see [ 12] for an overview), as well as its potentials to provide cheap electricity in large 
| quantities, available practically everywhere in the world |fT3l IT31 . 

Airborne wind energy generators exploit the motion of wings flying fast in the so-called crosswind conditions, i.e. roughly 
perpendicularly to the wind flow, and linked to the ground by flexible lines. The aerodynamic forces generated by the 
wing are then converted into electricity using one of several possible approaches 11121 . Control engineering plays a crucial 
role in airborne wind energy technologies, since, differently from conventional wind energy based on wind turbines, there 
is no passive, rigid structure that constrains the motion of the wing in the air and imposes its crosswind motion. Rather, 
this task has to be accomplished by an active control system that keeps track of the wing motion and issues suitable control 
inputs through actuators. In the last six years, the problem of control design for airborne wind energy generators has been 
studied by several research groups and companies, leading to a quite significant series of theoretical and numerical studies 
(see e.g. ifTTl l27l [8l [141 151 [151) as well as experimental tests, of which few works in the literature report measured data 

El El El. 

All of the mentioned control approaches rely on the availability of a series of variables to be used for feedback, most 
notably the wing's three-dimensional position and speed. The problem of estimating with sufficient accuracy and limited 
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lag these quantities is therefore of paramount importance in the field, however in the literature there is a lack of works 
on this topic, highlighting the specific issues that need to be addressed and providing either numerical or experimental 
results. 

In this paper, we contribute to fill this gap by providing a formulation of the estimation problem and by analyzing two 
different sensor setups and the related sensor fusion algorithms. In particular, the first sensor setup implies the use of 
a commercial Inertial Measurement Unit (IMU), installed on the wing and equipped with accelerometers, gyroscopes, 
magnetometers, GPS, and a barometer. We propose two algorithms for this sensor setup: the first approach is a quite 
standard linear observer, while the second one introduces a nonlinear correction to account for the specific kinematic 
constraints of the considered application. In the second sensor setup, we consider the use of a line angle measurement 
system in addition to the accelerometers, gyroscopes and magnetometers. The presented study is focused on the estimation 
of the wing's position and of its "speed angle", which is a variable well-suited for automatic control strategies, when 
the wing is flying fast in crosswind conditions. We present experimental results obtained with a small-scale prototype 
operating with a fixed line length of about 30 m and equipped with an IMU and a line angle measurement system, which 
we developed and built ad-hoc for this application. The equations employed in the design of the observer are based purely 
on kinematics, hence they are exact (i.e. they are not affected by uncertainty in the system's parameters) and independent 
from system's characteristics like the wing area, aerodynamic efficiency, mass, etc. Therefore, the presented experimental 
results are representative also of systems with larger wing size and different design, including systems with one and two 
tethers and/or rigid wings. 

2 System description, model equations and sensors setup 
2.1 System description 

The considered system is based on a small-scale prototype built at the University of California, Santa Barbara, shown 
in Fig. [T](see [4 | for a short movie clip). A flexible wing (inflatable power kites are used in the experiments) is linked 




Figure 1: Small-scale prototype for the control of tethered wings built at the University of California, Santa Barbara. 

by three lines to a ground unit (GU). In normal flight conditions, the wing's trajectory evolves downwind with respect 
to the GU. The GU is fixed to the ground and it is equipped with actuators, able to influence the wing's path: a human 
operator or an automatic control algorithm commands the actuators in order to obtain "figure-eight" crosswind paths, 
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which maximize the generated forces Ifl2l . This paper is devoted to the problem of estimating in real-time, in the described 
flying conditions, the wing's position as well as its "speed angle" of the wing. In the next section, we provide a formal 
definition of these variables and we introduce the related notation and model equations. 

The prototype operates with a fixed length of the lines of r = 30 m, hence it is not able to generate electricity through the 
synchronized reeling out of the three lines, as it is typical in ground-based airborne wind energy systems (| 12 1). Instead, 
the energy needed for the operation of the system is drawn from batteries installed on the GU. The use of a fixed, short 
lines' length (as compared with the intended operating conditions of this kind of systems, ranging from 100 to 1000 m of 
lines' length) does not impair the significance of the presented study, rather it increases it, for ar least four reasons: first, 
airborne wind energy systems operate with very low line unrolling speed as compared with the other components of the 
wing's speed vector, so the behavior of the wing does not change significantly between fixed line length and variable line 
length; second, the optimal operation of airborne wind energy generators is obtained with a constant line speed, and the 
settings considered here can be regarded to as a particular case of such operating conditions, i.e. with constant line speed 
equal to zero; third, the use of fixed lines' length yields the largest possible forces for given wind conditions, hence the 
highest accelerations and angular velocities, making the estimation problem more challenging; finally, the use of short 
lines implies that the wing's path is contained in a small portion of the aerial space, with consequent much more frequent 
changes of direction and inversions of motion, again increasing the difficulty of the filtering problem. 



2.2 Model equations 

Before stating the model equations used to design our filtering algorithms, we need to introduce a series of right-handed 
reference systems, as well as the transformations to convert a given vector from one reference to another. A first, inertial 
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Figure 2: (a) Coordinate system G and position of the wing p; (b) local coordinate system L = (Ljv, Le, Ld) and 
spherical coordinates 9, <j>; (c) wing speed angle 7; (d) coordinate systems G and NED and angle in- 
coordinate system G = (X, Y, Z) (Fig. |2]-(a)) is centered at the GU location, with the X axis parallel to the ground, 
contained in the symmetry plane of the GU and pointing downwind towards the wing. The Z axis is perpendicular to 
the ground pointing upwards, and the Y axis forms a right-handed system. The wing position vector expressed in the 
reference system G is denoted by pa = [px ,Py ,Pz] T & K 3 , where px, py and pz are the scalar components of pc 
along the axes X, Y and Z, respectively, and T stands for the matrix transpose operation. The wing's position can be 
also expressed in the spherical coordinates 9, <fi (Fig. |2}(b)), where 9 £ [0, ?] is the angle between the (X, Y) plane and 
vector pc and <f> e [— tt, tt] is the angle between the X axis and the projection of pc onto the (X, Y) plane, taken to be 
positive for a positive rotation around the Z axis. In particular, we have: 
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where r is the distance between the GU and the wing, equal to the lines' length in our setup. We note that the angle 9 
considered here is the complement to ^ with respect to the angle employed in previous works (see e.g. IS)), and that the 
9, 4> angles are different from those used in [1 lj . 

Angles 9, <f> also define a second, non-inertial Cartesian coordinate system, L = (Ljy, Le, Ld), which we call "local". 
Vector pa (as well as any other vector) can be expressed in system L by using the rotation matrix R G ^ L : 

PL = Rg^lPG, (2) 



R G - 



— sin (9) cos ((/>) — sin (9) sin ((f)) cos (9) 

— sin {(f)) cos ((f)) 

— cos ((f)) cos (9) — sin ((f>) cos (9) — sin (9) 



(3) 



and the inverse transformation is given by the rotation matrix R, _._ = R„ 1 r = RfZ , . The axes Ln and Le define the 
tangent plane at point pa to the sphere of radius r, on which the wing's trajectory is confined, and they can be interpreted 
as local north and east direction relative to such a sphere. Hence, Lb is always parallel to the ground (i.e. to the (X, Y) 
plane). The axis Le thus represents the local down, pointing from the wing to the center of the sphere (i.e. to the GU 
location). 

We further consider the wing's coordinate system, K = (K x , K y , K z ), and standard geographical North East Down 
system, NED = (N,E,D). System K is centered at pq, it is non-inertial and fixed with respect to the wing, i.e. 
it provides the wing's orientation. The rotation matrix R NED ^ G is used to express a vector in 7VLZ3-coordinates into 
G-coordinates: 

/cosOg) smOc) \ 
Rnbd^g= sin (<Pg) -cosOg) , (4) 
V 0-1/ 

where (f)Q e [0, 2ir] is the angle between the N axis and the X axis, measured by a positive rotation around the D axis 
(see Fig. EKd)). Moreover, a vector expressed in the K system can be converted into the NED system by means of the 
rotation matrix R K ^ NED (<?), where q — [qi, q 2 , qz, qi\ T £ K 4 is the quaternion defining the relative orientation between 
K and NED (see e.g. USUI): 



(5) 



R'K—rN ED \Q) 

'%i %293 - 9194) 2(9294 + qmT 

2(^293 + 9i94) 2(gf + qf)-l 2(q 3 q 4 -q 1 q 2 ) 
K 2(q 2 q 4 - qiq 3 ) 2{q 3 q 4 + qi q 2 ) 2(gf + qf) - l J 

The matrix R K ^ G = R NED ^ G R k ^ned can used to express a vector in K coordinates into G coordinates. 
Finally, we define the wing's speed angle 7 e [— it, ir] as follows (see Fig. |2-(c)): 

7 = arctan 2 (pl e ,Pl n ), (6) 

where arctari2 (pl e ,Pl n ) € [— 7r, it] is the 4-quadrant arc tangent function and pl e , Pl n are proportional to the sine 
and cosine components of 7, respectively. The angle 7 provides the angle between the wing's speed vector and the local 
north axis, Ln, measured positive for a positive rotation around the local down Ld. We note that, since the considered 
system has fixed line length, vectorp is always contained in the (Ln, Le) plane. In case of variable line length, equation 
(O has to be adjusted by considering the projection of p onto (Ljy, Le)- Angle 7 is particularly important as feedback 
variable for automatic control algorithms, since it can be easily linked to the wing's path: as an example, 7 = 0, 7 = \ 
and 7 = 7r indicate that the wing is moving, respectively, towards the local north, parallel to the ground towards the local 
east, or towards the ground. 

We are now in position to introduce the model equations we consider in this work. These equations are based entirely 
on the system's kinematics, i.e. they are basically obtained by differentiating twice the vector pc(t) with respect to the 
time variable t. Kinematic equations bring two important advantages: 1) they provide an exact model, i.e. there are no 
neglected dynamics even for flexible wings, and 2) they do not depend on any of the system's characteristics, like mass, 
shape, moments of inertia and aerodynamics of the wing. Indeed these features influence the motion of the wing through 
complex, infinite dimensional nonlinear dynamics, whose inputs are the steering command given by the GU and the 
(unmeasured) wind, yet such dynamics are totally irrelevant for our scope if a measure or estimate of pc(t) is available. 
However, while the kinematics in the inertial frame G are given by linear operators (i.e. derivatives), we still end up with 
nonlinear model equations, due to the fact that we measure the involved variables in different reference frames. In order 
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to facilitate the observer design, we split the model into five interconnected subsystems, as shown in Fig. [3] The involved 
variables are the wing position, speed and acceleration in the reference G, Pcj(t), Pa(t), Pg(*)> the wing acceleration 
in the reference K, PK(t), the quaternion q(t), finally the wing angular velocity vector in the reference K, cDfj-(t). The 
model equations are: 
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where O3 is a 3 x 3 matrix of zeros, and J3 is the 3x3 identity matrix. Thus, the considered model has px , ^>k (t) as inputs 



9,4, 



Figure 3: Scheme of the considered kinematic model. 

and pc(t), Pc(t), q(t) as states, and it is composed by a linear dynamical system, C, a nonlinear one, TV, and three static 
nonlinear functions, fx, fa, f%. The only involved parameters are the constant angle cj>Q, giving the orientation of the GU 
with respect to the geographical north, and the length r of the lines. Both these parameters can be measured accurately 
and are assumed to be known exactly here. In the next section, we will describe the different considered sensor setups and 
indicate the corresponding measured variables. 



2.3 Sensors setup 

Several sensors are installed on the GU and onboard the wing, giving a range of different possibilities for the design of 
observers forpG(i), PG{t), and j(t). A GPS and three magnetometers are installed on the GU, providing its geographical 
location in NED as well as the angle <j)Q. Moreover, the GU is equipped with a line angle measurement system (shown 
in Fig. 0J, which provides a direct measurement of the angle of the main line connecting the wing to the GU. This sensor 
has been developed and built ad-hoc for this application and it employs two incremental encoders. The encoders are 
connected to the center line via a metal rod and a small pulley which allow the center line to slide with low friction with 
respect to the rod while the wing is flying. The first encoder measures the angle between the metal rod and the (X, Y) 
plane, while the second one measures the angle between the X axis and the projection of the metal rod onto (X, Y). Both 
encoders have a resolution of la( \^ j e 400 counts for each full revolution. Driven by the movement of the center line, 
the metal rod rotates following the wing's path, and the encoders measure such rotations in terms of the above mentioned 
angles. More specifically, the readings from the two incremental encoders, denoted as 9b and 4>b, can be converted into 
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Figure 4: Line angle sensor 



measurements of 9 and via the following equations (see Fig. |5): 



L=^L{ + Li, (8a) 
9b = @b — arctan2 ( — ) (8b) 

l[= Lain (9 B ) (8c) 
l 2 = L cos (9 B ) cos (0b) — ?2 (8d) 
2 3 = L cos sin (c^b) (8e) 

6» = arctan ( \ l + i | (8f) 

= arctan 2 (l' 3 , / 2 ), (8g) 

where Li and L2 are, respectively, the lengths of the metal rod and the pulley on top of it, and £1, I2 are the distances 
between the encoders' position P2 and the attachment point of the center line on the GU, Pi . All quantities L\, L2, h, h 
are fixed and known. 

As regards onboard sensors, the wing is equipped with an Inertial Measurement Unit (IMU) manufactured by SBG 
Systems®, comprising three accelerometers, three gyroscopes, three magnetometers, a barometer, and a GPS. A 900 
MHz radio is also installed on the wing and transmits the measurements collected by the IMU to a receiver on the GU, at 
a rate of 1 15200 Baud. The accelerometers have a measurement range of ±5 g (where g is the gravity constant) and noise 
densityof2.5xl0- 4 g/%/Hz. The gyroscopes have a measurement range of ±300 °/s and noise density of 0.05 ° j (s\/Hz). 
The magnetometers have a measurement range of ±1.2 Gauss and noise density of 1 10~ 5 x Gauss/ \/Hz. The barometer 
provides a measurement of barometric altitude with a resolution of 20 cm. All sensors mentioned so far provide a new 
measurement with a sampling time of T s = 0.02 s. Finally, the onboard GPS has a nominal horizontal accuracy (i.e. on 
the (N, E) plane) of ±2.5 m and a sampling time of 0.25 s. 
In the following, we will implicitly assume that: 

• the line angle sensor provides a measure of the 9, <fi angles; 

• the position measured by the onboard GPS coincides with pc\ 

• the barometer provides a measure of pz; 

• the IMU is fixed with respect to the kite, hence the accelerometers, gyroscopes and magnetometers provide mea- 
surement of, respectively, fix, ^k, and q. 
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Figure 5: Diagram of the line angle sensor: (a) side view and (b) top view. Dash-dotted line: center line of the wing; thick 
solid lines: metal rod and pulley linking the encoders to the center line. 



These assumptions imply that the wing's lines are parallel to its position vector pq, that the position of the IMU coincides 
with pa and that the IMU itself is fixed to the wing, so that the measured accelerations and angular speeds corresponds to 
those of the wing. Given, the short length of the wing's lines, such hypotheses are reasonable in the considered application. 



3 Sensor Fusion for Airborne Wind Energy 

In this section we present different algorithms to estimate the wing's position and speed angle, pc(t) and j(t). In 
the following, we use the notation to indicate the noise-corrupted measurement of a given variable, e.g. pa is the 
measurement of the position vector seen from system G, pa- Moreover, we indicate with '"' the estimated variables, i.e. 
the output of the observer. 

As a first step, we exploit the structure of the considered kinematic model to separate the problem of estimating the 
wing orientation (i.e. the quaternion q(t)) from the one of estimating its position and speed. In fact, since the nonlinear 
system Af is known and the IMU provides a measurement of both its input and its state, i.e. unit) and q(t), an extended 
Kalman filter (EKF) approach can be adopted to obtain the value of q(t). The use of EKFs for the absolute orientation 
of rigid bodies using gyroscopes and magnetometers has been already studied in the literature (see e.g. |22j), and the 
IMU employed in this study already includes its own EKF for this purpose. We indicate such a filter as Nekf- With 
the estimate q(t) provided by Nekf, we can then compute the rotation matrix RK^NED{q{t)) by using (0. Then, we 
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compute an estimate po(t) of pc(t) through the static nonlinear function /i given by: 

f 1 :p G (t) = R 

NED—+G K N E D 



(9) 



/i is obtained from fi (TTct by using the estimated rotation matrix RK^NED{q(t)) instead of the "true" one, and then by 
removing the bias given by the gravity acceleration g. The latter is present in our setup due to the characteristics of the 
employed MEMS accelero meters. The approach described so far is used as a preliminary step for any of the observers 
described in the following sections, whose goal is to estimate the wing's position and speed. In particular, consider that 

the estimate pG(t) can be re-written as: 

p G (t)=PG(t)+rj p (t), (10) 

where fj p (t) g R 3 is some estimation error that we consider as process noise. Then, the estimate pc(t) can be regarded to 
as a noise-corrupted input for the linear system C (see Fig. |3). Following this idea, we will use a steady-state Kalman filter, 

based on C, to compute pa(k) and pc(k)- In order to do so, we discretize C by forward difference with the employed 
sampling time T s : 



p G (k+l) 
Pa(k) = [I 



h 
0, 



T a I 3 
h 



P_c(k) 
Po(k) 




pa{k) + ff p (k) 



3 



A 

Pa{k) 
Pa(k) 



(11) 



+ Vm(k), 



where k g Z is the discrete time variable and ff m E M. 3 represents the measurement noise. In classical Kalman filtering 
theory, ff p and rj m are assumed to be independent white Gaussian processes with covariance matrices Q and R, respec- 
tively |26l . Here, the measurement and process noises can be reasonably assumed to be independent, since they pertain 

to completely different sensors (accelerometers, gyroscopes and magnetometers for the input pc{k) and either GPS and 
barometer or line angle sensor for the output pc{k)), yet they are not white Gaussian processes, since the measurements 
are the result of algorithms like Nek f or the GPS algorithm. However, the matrices Q and R can still be used as tuning 
knobs in the design of the linear observer: a larger Q with respect to R will yield an observer that relies more on the 
measured output and carries out more aggressive corrections of the estimates, viceversa a larger R will yield a smoother 
course of the estimated variables, with higher discrepancy between the measured and estimated outputs. 
The steady-state discrete Kalman Filter is then obtained by the following equations 
Time update ("Predict"): 



Vpik) 



Measurement Update ("Correct"): 



Pc(k) 
Pc(k) 



= A 



Pp{k) 



Pp(k - I) 



Bp G {k) 



K p G (k) - C 



Pp( k ) 
Pc(k) 



(12) 



(13) 



where p G (k), p G (k) are the a-priori state estimates at step k given the knowledge of the process prior to step k, while 

Pa(k), PG{k) are the a-posteriori state estimates at step k given the knowledge of the measurement pcik). The steady- 
state Kalman gain K is computed off-line as: 

K = AP^iCP^C 7 + R)- 1 , (14) 
where P^ satisfies the following Algebraic Riccati Equation: 

Poo = AP^A 7 - AP OQ C T {CP QO C T + R^CPooA 1, + Q. (15) 
We will now describe the three different approaches considered in this study. All three of them employ the described 
procedure to obtain pc(k) and then the Kalman filter (fT2l-([T3l to estimate position and speed, but they use different 
strategies and sensors to obtain the position measurement pc(fc). A scheme of each of the three approaches is shown in 
Fig. [6] where we indicate the Kalman filter as Ckf- 
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Figure 6: Schemes of the designed observers for the first, second and third approach ((a), (b) and (c), respectively). 



3.1 First approach: GPS and barometer 

We design a first observer (Fig. |6}(a)) by using the GPS to obtain the measurements px, Py, and the barometer for pz- 
Hence, the two sensors together provide the position measurement pc(k) to be used with the Kalman filter (fT2l>-(ll~3l>. 
While the GPS can also provide a measurement of the altitude, the related error (about ±50 m in our setup) is too large 
to be used in the considered application. The barometer, on the other hand, is quite accurate provided that an initial 
tuning procedure is carried out for each test to remove the bias induced by the changing weather conditions. The GPS 
measurements are obtained by computing the difference of the readings between the onboard GPS and the one placed on 
the GU. Since the accelerometers have a sampling frequency of 50 Hz, while the GPS provides measurements at a 4 Hz 
rate, we use a multi-rate Kalman filter in this first approach, which consists in performing the time update (fT2] i at each 
sampling time of T s = 0.02 s, and the measurement update dl~3l every 0.25 s, i.e. when the new GPS measurement is 
available. 



3.2 Second approach: GPS and barometer with geometric correction 

The filter designed in the second approach employs the same sensors as the previous one, but it also carries out a correction 
of the measured position on the (X, Y) plane given by the GPS. Such a correction is based on two observations: 1) the 
wing motion is constrained on a sphere of known radius r, 2) the measurements px , Py given by the GPS are less reliable 
than the measurement pz provided by the barometer. Using Z, by inverting (l7dl i we compute a measure 9 of angle 9, 
then, we correct px, Py by re-scaling them in order to match with the line length r projected onto (X, Y), i.e. r cos(#) 
(compare ([TJ and Fig. |2}(b)): 

rcos(^) 
PX,c = VX—f== 

VPx+Py (If.) 
~ rcos(6>) K UJ 

Then, we use the vector pc, c = [Px,c, Py,cPz] T as position measurement for the multi-rate Kalman filter. This second 
approach is shown in Fig. |6r(b), where the correction ( fT6b is indicated as Me- 
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3.3 Third approach: Line angle sensor 

The third observer we consider (Fig. |6}(c)) employs the line angle sensor to obtain the position measurement. Hence, in 
this approach we move away from using the GPS and barometer measurements. By inverting (l7dl > we have: 

Px(k) = r cos(0(fc)) cos(0(fc)) 

p Y (k) = rcos(#(fc))sin(>(fc)) (17) 
p z (k) = rsm(9{k)), 

where 0(k), 4>(k) are the measurements provided at 50 Hz by the line angle sensor. Due to the use of incremental en- 
coders, the measurement noise on 9(k), 4>{k) is basically given just by the quantization error resulting from the encoders' 
resolution, i.e. rad. 



3.4 7 filter 

In order to obtain an estimate of the speed angle 7 ©, we first use dTdb to obtain the angles 9, <j> from the filtered position 
Pa- The latter is obtained, together with the speed estimate pQ, using one of the observers presented in the previous 



sections. Then, we estimate the rotation matrix R a _ il 
in the L frame, p^, as: 

PL 



by using (O, and we compute an estimate of the speed vector 



R G ^ L {6, <j>)p G . 



(18) 



Using pl, we compute a first "measurement" 7 ~ 7 by simply applying its definition @. However, 7 results to be 
affected by high-frequency oscillations (see section [4] for an example with experimental data), hence it is not suitable to 
be used for feedback control. In order to have a smooth estimate, we design another observer based on the following 
state-space model: 



(19) 
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In particular, we design a standard Luenberger observer ifTOl to obtain the estimates 7 and 7. 

Clearly, the three different approaches to estimate the wing's position and speed yield estimates of 7 with different 
accuracy. In the next section, we discuss the performance obtained with the considered filters in experimental tests. 



4 Experimental results 

We designed three different observers for pa, according to the approaches described in section[3] and consequently three 
observers for 7. In this section, we use the superscripts 1 , 2 , 3 to indicate the filtered position obtained with the approach 
of section [3~Tll3.2l and l3.31 respectively. For the first and second approach, we used the following covariance matrices: 
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while for the third approach we used: 
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The gain for the Luenberger observer used to estimate 7 is equal to K~ f = [0.06, 0.6] T . The length of the lines is r = 30 m. 
We employed different wings in the experiments, in order to demonstrate that the proposed approach is independent from 
system's features like wing area, mass, efficiency, etc. In particular, we used three different Airush One® inflatable power 
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Figure 7: Experimental results. Estimates of the position along the X axis, px{t), obtained with the first (dotted line), 
second (dashed) and third (solid) approaches. 
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Figure 8: Experimental results. Estimates of the position along the Z axis, pz(t), obtained with either the first and second 
(dashed line) or the third (solid) approach. 

kites, with 12 m 2 , 9 m 2 and 6 m 2 area. 

The third approach results to be the most accurate one, thanks to the use of the line angle sensor, which practically yields 
exact position readings. Hence, in this section we take the estimates and 7 3 , obtained with the third approach, as 
"true values" and we evaluate the deviations resulting from the first two approaches with respect to the third one. An 
example of the courses of the position component along the X axis, px, is shown in Fig. [7] It can be noted that the 
error between p x , p x and p x can be quite high (of the order of 10 m), with p 2 x being generally better than p x . This 
shows that the correction ( TToT ) we introduced in the second approach indeed improves the estimate with respect to using 
the plain GPS reading. Moreover, both p x and p\ are affected by some delay with respect to p x . The estimates of 
vertical position result to be more accurate thanks to the use of the barometric altitude instead of the GPS, see Fig. [8]for 
an example (we note that for the vertical position pz, the first two approaches are equivalent hence resulting in the same 
filtered values). Table Q] shows the root mean square errors (RMSE) between the first two approaches and the third one, 
for all of the estimated variables. The results of Table Q] account for about 5 hours of operation at different wind speeds. 
Considering that a single figure-eight loop has a period of 5-6 seconds, the data in the table correspond to about 3300 
complete figure-eight paths. The results in the table confirm that the second approach is generally better than the first one 
in estimating the wing's position, especially at larger wind speeds (see TableQ]for 3.5 m/s and above), while the two give 
similar errors on the 7 estimates, which deviate significantly from the third approach. Moreover, it can be noted that the 
accuracy of the first two approaches get worse as the wind speed is higher. Since the wing's speed depends linearly on 
the wind speed (see e.g. 1 14 1), this results indicate that the GPS performance get worse as the wing movements get faster 
and changes of direction are more frequent. Indeed the provided GPS accuracy (±2.5 m in the (TAT, E) plane) is valid for 
steady state conditions only. 
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Figure 9: Experimental results. Courses of j(t) (dashed line) and j(t) obtained with the third approach. Wing area: 12 
m 2 . 
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Figure 10: Experimental results. Path of the wing during one figure-eight cycle in the (<f>, 9) plane (dotted lines) and 
gradients obtained with the estimated speed angle j(t) obtained with the third approach (arrows). Wing area: 12 m 2 . 

Fig. |9]shows, relative to the third approach, the smoothing effect obtained with the Luenberger observer for the j 3 angle 
with respect to the values of 7 3 . A typical example of figure-eight path in the (0, 9) plane, obtained during the tests, is 

shown in Fig. [TO] together with the gradient vectors [(f), 9] T estimated using the values of 7 3 : it can be noted that the 
estimated speed angle is quite accurate, with the gradient being almost always tangent to the wing's flying path. The 
length of the arrows in Fig. [TU]is proportional to the wing's speed with respect to the ground. Similar good results have 
been obtained with all the three different wings, and in different wind conditions (see Fig. QT|for some examples), hence 
showing the robustness of the approach deriving from the use of kinematic equations to design the filters. These features 
make the use of the proposed approach and of 7 3 most suited for feedback control, see [4| for a short movie clip concerned 
with autonomous flights carried out using such estimation algorithm. 



Table 1: Experimental results. Root mean square errors between the estimates obtained with either the first or second 
approach and the third one. 
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pi 2 nd 

5.61 2.36 
7.28 4.16 
0.74 0.74 
1.70 0.80 


,si 2 " J 

7.74 5.07 
10.32 10.16 
1.71 1.71 
1.35 1.44 
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10.67 6.64 
14.18 12.56 
2.13 2.13 
1.52 1.71 


pi 2 L,d 

12.81 6.96 
12.77 11.16 
2.75 2.75 
1.50 1.61 


pi 2 nd 

13.30 7.80 
13.27 10.71 
3.47 3.47 
1.54 1.72 


pi 2 nd 

17.18 9.71 
14.97 10.98 
4.58 4.58 
1.71 1.97 


pi 2 nd 

28.91 12.45 
18.02 11.12 
6.36 6.36 
1.90 2.32 


pi 2 nd 

47.21 14.50 
16.47 10.31 
8.44 8.44 
2.40 2.37 


pi 2 " d 

55.39 14.55 
13.35 9.45 
9.30 9.30 
2.67 2.24 


pi 2 nd 

56.61 14.93 
12.83 9.61 
9.15 9.15 
2.65 2.35 


pi 2» d 

56.39 15.23 
13.46 9.81 
9.32 9.32 
2.96 2.39 
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5 Conclusions 



We presented a study on the problem of sensor fusion for airborne wind energy systems. We designed three algorithms 
to estimate the wing's position and speed angle, using different sensors, and we applied them in experimental tests. The 
presented results allowed to draw some conclusions on the efficacy of the different employed sensors for this kind of 
application: for relatively short line length, GPS sensors are not suited to track the highly dynamical movement and 
frequent change of direction of the wing, while line angle sensors results to be very effective. This situation is likely to 
change for longer line length and larger paths, when the wing keeps a constant heading for longer time: in this context, 
the GPS performance should improve, while the line angle measurement might be less accurate due to eventual line sag. 
More expensive GPS sensors and differential GPS systems should also be able to obtain better performance, although 
for short line length the accuracy and bandwidth of a line angle sensor cannot be matched. The presented algorithms are 
based on a kinematic model of the system, which is exact and independent from the features of the particular employed 
wing. Hence, the obtained results are significant also for other wing sizes and for rigid wings. 
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